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The feasibility of a mean-field dynamo in nonhelical turbulence with superimposed linear shear 
is studied numerically in elongated shearing boxes. Exponential growth of magnetic field at scales 
much larger than the outer scale of the turbulence is found. The charateristic scale of the field is 
l-g oc S' -1 '' 2 and growth rate is 7 oc S, where S is the shearing rate. This newly discovered shear 
dynamo effect potentially represents a very generic mechanism for generating large-scale magnetic 
fields in a broad class of astrophysical systems with spatially coherent mean flows. 
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Introduction. Understanding the origin of cosmic 
magnetism is one of the fundamental theoretical chal- 
lenges in astrophysics. The turbulent motions of the 
plasmas that make up most astrophysical objects are be- 
lieved to be responsible for the generation of the mag- 
netic field. In particular, a generic property of the tur- 
bulence of conducting fluid is to amplify exponentially 
magnetic fluctuations at the turbulence scales or smaller 
via the fluctuation dynamo effect [H, 0, S|- A distinct 
problem is to explain the observed presence in most as- 
trophysical bodies of magnetic fields spatially coherent 
at scales larger than the outer scale of the turbulence 
(mean fields). Nonhelical homogeneous isotropic turbu- 
lence on its own cannot give rise to a mean field. What 
are then the large-scale properties that must be present 
in a turbulent system for such a field to be generated? 
Mean-field dynamo theories [3] have identified a number 
of amplification mechanisms. We know, e.g., that non- 
zero net helicity (often combined with rotation in real 
systems) is sufficient to produce mean fields, but is it 
necessary? 

Perhaps the most common large-scale feature is a mean 
velocity shear. It is present, e.g., in stellar interiors 
accretion disks 6] , galaxies (in particular, irregular ones 
0), and liquid- metal laboratory dynamos all of which 
host both large-scale (mean) and small-scale (fluctuat- 
ing) magnetic fields. A number of theories have proposed 
that a mere combination of turbulence and shear could 
give rise to a mean-field dynamo: e.g. , the shear-current 
effect [ij], the stochastic a effect shear amplifica- 

tion of small-scale-dynamo-generated field ll| , negative- 
diffusivity-type theories [lij . Ref. @, which used the 
r-approximation closure, provoked a particular debate 
because it seemed to contradict the rigorous mean field 
theory based on the second-order correlation approxima- 
tion (SOCA), which ruled out the shear dynamo [HI]. 



However, the SOCA is only valid in the limit either 
of low hydrodynamic and magnetic Reynolds numbers, 
Re,Rm <C 1, or short velocity correlation times The 
real turbulent systems are in neither of these limits, and 
the hope that some of the results qualitatively carry over 
has had to be backed up by numerical evidence Il4jl and 
by intuitive physical field-amplification scenarios [la ] . In 
the absence of a compelling physical argument for or 
against the shear dynamo or of a rigorous method for 
proceeding analytically, a numerical experiment is over- 
due. Here we report the first such numerical experiment, 
which supports the existence of the shear dynamo. 

Numerical Set Up. We consider the incompressible 
magnetohydrodynamics (MHD) with a background lin- 
ear shear flow U = —Sxy and a white- noise nonhelical 
random homogeneous isotropic body force f: 
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where u and B are the velocity and magnetic fields, 
dj dt = dt — Sxdy + u • V, the density p ~ 1 , and the pres- 
sure p is determined by the incompressibility condition 
V • u = 0. These equations are solved with shear-periodic 
boundary conditions by a Lagrangian spectral method 
[l6| . When the imposed shear S is weak compared to the 
turnover rate of the turbulent motions, the growth of the 
mean (large-scale) field can only be detected if the size of 
the computational domain is much larger than the turbu- 
lence scale Iq- In general, this, together with the necessity 
to run the simulations for very long times, requires un- 
affordable amounts of computing power. We circumvent 
this problem by using computational boxes with large as- 
pect ratios, L x x L y x L Zl where L z L x = L y . The 
units of length and time are fixed by setting L x = L y = 1 
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TABLE I: Index of runs 
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and the mean forcing power e = (u • f ) = 1 (this can 
be controlled because the forcing is white- noise) . The 
average forcing scale is Iq = 1/3, i.e., the energy is in- 
jected randomly in the wave-number shell centered at 
ko/2ir = 3. The resulting root-mean-square velocity field 
is u lms = (it 2 ) 1 / 2 ~ 1 ± 0.3, so the typical turnover 
rate of the turbulent motions is M rm s/^o ~ 3. We study 
five values of the shear S = 2,1, 1/2, 1/4, 1/8 < u lms /lo 
(weak shear is used in order for the mean and fluctuat- 
ing fields to be clearly distinguishable from each other 
via small-scale averaging; see Eq. ([3]) below). The vis- 
cosity and magnetic diffusivity are v = r\ = 1CP 2 , so 
Rm = Re = u rms /kois ~ 5. The resolution requirements 
are consequently not large: it suffices to have 32 x 32 col- 
location points in the (x, y) plane. In the z direction, we 
use resolutions between 256 and 4096 collocation points 
for L z = 8, . . . , 128, depending on S (Tab. HJ. 

Strictly speaking, we cannot speak about turbulence 
with such low Re. However, a developed inertial range 
is not important for mean field dynamos: it is sufficient 
that a stochastic velocity field with Re > 1 is present 21 1 . 
In our simulations, Rm is subcritical with respect to the 
fluctuation dynamo @, 0| , so any field growth we detect 
is due purely to a mean-field dynamo. Note, however, 
that since Rm > 1, turbulent tangling of the mean field 
generates small-scale magnetic fluctuations whose energy 
is in general not smaller than that of the mean field [3] . 

Results. All runs are initialized with a random, zero- 
mean, dynamically weak ((B 2 ) = 10~ 20 ) magnetic field. 
The field grows exponentially with time at all values of 5 
studied, provided the computational box is sufficiently 
long. For each S, we consider the growth rate 7 of 
B rms = (i? 2 ) 1 / 2 to be converged if it stays approximately 
the same when L z is doubled (Fig.[T]). That we are able 
to find such values means that 7 is asymptotically inde- 
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FIG. 1: Evolution of u lms (upper panel) and -B rms (lower 
panel) for 5 = 1 and four values L z — 8, 16, 32, 64. 

pendent of L z (the dependence on L x and L y should also 
be studied but that is currently too expensive computa- 
tionally). The field eventually grows to a dynamically 
strong saturated level. Here we concentrate on the kine- 
matic (weak-field) regime and leave the properties of the 
saturated state to a future study. 

Fig. [5] shows that, in the range of shears studied, the 
growth rate of -B r ms appears to increase linearly with S, 
7«S. This is a somewhat unexpected result from the 
theoretical point of view because the shear-current ef- 
fect [9( , as well as most other mean- field theories quoted 
above predict 7 cx S 2 for the fastest-growing mode. We 
cannot, however, exclude the possibility that the S 2 scal- 
ing may be asymptotically recovered at even smaller S. 

That the growing field is large-scale is obvious already 
from the visualization of the field: the large-scale as- 
dependent modulation is evident against the turbulence- 
scale structure (Fig. [3]). We isolate this large-scale de- 
pendence on z by low-pass filtering in Fourier space: 



B(z)= B(fc, = 0,fc y = 0,fc z )e 

|fcz/27T|<l 



(3) 



Note that since V • B = 0, B z = 0. This procedure av- 
erages out the small-scale structure and brings out the 
growing large-scale field in a clear way (Fig. [3]). In all 
cases, the root-mean-square values of B x and B y grow 
exponentially with the same rate 7 as B lms . We find 
\B y \ > \B X \, which is expected because the shear sys- 
tematically converts B x into B y [Eq. @]. The ratio 
\B y /B x \ S/^f is approximately constant in time and 
its average is independent of S (Tab. HJ, which is consis- 
tent with 7 cx S established above. 

Examining Fig.[3j we see that the magnetic field grows 
in large random patches. In time, they move around (in 
z) and change shape in an apparently random fashion. 
This means that the mean field is not strictly speaking 
an "eigenmode" with a spatial profile constant in time, 
although we found the correlation time scale of its evo- 
lution to be much longer than either the turnover time 
(^o/i*rms) or the shear time (S* _1 ). In order to describe 
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FIG. 2: Growth rates 7 of B lms for all runs (Tab. HJ). The 
dotted line shows the slope corresponding to 7 oc 5*. 

the spatial structure of the mean field in a systematic 
way, we define the time-averaged characteristic scale l^r: 
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Here and in all other cases, the time average is taken over 
the exponential-growth (kinematic) period (ii,<2) of the 
field evolution. The derivatives are calculated in Fourier 
space. The values of % are given in Tab. U and plotted 
vs. S in Fig. |4j As the shear is decreased, hg increases 
and is well matched by the scaling % oc S~ x l 2 . 

This scaling is again at odds with the mean-ficld-theory 
prediction l-g oc S" -1 9], but the following simple argu- 
ment shows that it is consistent with 7 oc S and suggests 
the possible form that a mean-held theory of the dynamo 
reported here may take. Let us conjecture the mean-field 
equations in the following (standard) model form 



d t B x 
d t B v 



-Vrk 2 z B x 
-VTk 2 z B y 



ABy, 

SB X , 



(5) 
(6) 



where rjx ~ Wrms^o is the turbulent diffusivity and A is 
some operator that closes the dynamo loop, (the chal- 
lenge of mean-field theory is to find A). The growth 



FIG. 4: The characteristic scale of the magnetic field [Eq. Q] 
for all runs. The dotted line showes the slope S _1//2 . 

rate is 7 = V SA — rjxk 2 . For the fastest-growing mode, 
the two terms in this expression are comparable, so, if 
k z ~ lz± ~ ^(SVurms) 172 , then A ~ S and 7 - S 
(in contrast, some mean- field theories predict A ~ Sl\k 2 z 
[iEi, so k z oc S and 7 ex S 2 ). 

Finally, in Fig. [5l we show the spectrum of magnetic 
energy during the growth stage. It is strongly peaked at 
large scales (k z lo -C 1), but also shows that the mean 
field is tangled by the turbulence to produce a significant 
amount of magnetic energy at the turbulence scales. The 
presence of this tangled component is likely to be essen- 
tial in the shear dynamo effect in that the small-scale 
field could be continuously resupplying the field B x for 
amplification by shear The key to a successful the- 
ory of the shear dynamo may be to find the way in which 
the symmetry is broken to make this resupply system- 
atically favorable to the exponential growth of the mean 
field (as, e.g., implicitly attempted in Ref. Q). 

Effect of Shear on Velocity Field. It has been sug- 
gested that a mechanism similar to the shear dynamo 
may also produce large-scale velocity structures ("vor- 
ticity dynamo" [l7|). The velocity does indeed develop 
large fluctuations that are energetically comparable to 
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FIG. 3: Snapshots of u y (upper panel) and B y (lower panel) 
taken in an (y, z) cross-section of the L z = 16 run for S = 1 
(the forcing scale is lo = 1/3). Underneath the snapshots are 
plots of u y (z), u x (z) (upper panel) and B y (z), B x (z) (lower 
panel). Here u(z) is defined similarly to B(z) [Eq. J3)]. 



■3 



10" 



10- 



10" 



10" 



10" 





1 


1 ; ; 

^i::: ! 1 






\ — "S- '*. ~ 

, . !% n 


; L z 

■ ■— L z 


= 8 
= 16 
= 32 


rms forcing , \>Y- 
wave number ! ViV* 


\ — 


(fc' ms /27r = y/3) i \$S 1 

V. m 


: — L z 


= 64 
1 



0.01 0.1 1 10 

k z /2ir 

FIG. 5: Normalized one-dimensional spectra of the magnetic 
energy, M(k z ) = Ylh x ,h l B ( fc *i k z )\ 2 /(B 2 ) (time averaged 
over the growth stage) for S = 1 and L z = 8, 16, 32, 64. The 
four graphs demonstrate that, as L z is increased, a large-scale 
spatial structure independent of the box length emerges. 
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the small-scale turbulence, last for long times (Fig. [J) 
and are spatially coherent on scales similar to those of 
the magnetic field (Fig. [2]). The large-scale structure 
forms mainly in u y (corresponding to large-scale vorticity, 
ZJ X = —duy/dz). As the forcing always seeds some large- 
scale vorticity that is not infinitcsimally small, the "vor- 
ticity dynamo" always operates in the nonlinear regime. 
Its detailed study is outside the scope of this Letter. 

The growth of the magnetic field does not seem to 
be strongly correlated with the evolution of the shear- 
generated large-scale velocity structures (compare, e.g., 
the time evolution of u [ms and B rms in Fig. [1} . 

It is well known that the presence of shear can lead 
to nonlinear destabilization of finite velocity fluctuations 
and formation of shear-driven turbulence whose outer 
scale is the scale of the shear (in simulations with a lin- 
ear shear, the box scale). This does indeed happen in 
our simulations when S is too strong or the box is too 
long. The quantitative signature of this regime is that 
the power input from the shear in Eq. |T]), (u x u y )S, ex- 
ceeds the forcing power e = (u • f ) . Wc avoid this regime 
to isolate the mean-field generation effect, which requires 
a scale separation between the turbulence and the mean 
field. In all runs reported here, \(u x u y )S\ <C e. We note 
that the large upward fluctuations of u rms [Fig. [T] are 
not accompanied by a significant change in (u x u y )S, so 
the large-scale velocity structures appear to feed on the 
forcing power, not on the power extracted from the shear. 

Discussion. We have found that a large-scale mag- 
netic field grows exponentially in long sheared boxes with 
forced small-scale nonhelical turbulence. In the param- 
eter range we have studied, the growth rate appears to 
scale as 7 oc S 1 , the spatial scale of the field l-g cx S^ 1 ^ 2 
and B y /B x ~ const > 1. These properties do not seem 
to fit any of the existing theories. Our results do, how- 
ever, lend credence to the concept of a shear dynamo and 
thus provide motivation for further theoretical effort. 

To our knowledge, this is the first demonstration of the 
shear dynamo effect in a dedicated numerical experiment. 
In an earlier unpublished study we obtained similar re- 
sults using PENCIL code (a compressible finite-difference 
code in contrast to the spectral one used above), so the 
amplification effect appears to be numerically robust. We 
note that there have been earlier indications of nonheli- 
cal turbulence amplifying large-scale magnetic field in the 
presence of a large-scale shear associated with mean flows 
in numerical experiments t hat used constant-in-time si- 
nusoidal forcing functions 
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Another example of 
large-scale magnetic fields generated by a combination of 
nonhelical turbulence and a mean flow is the numerical 
experiments with Taylor-Green forcing [20| . One might 
speculate that the shear provided by the mean flow in 
such systems could act in a way qualitatively similar to 
a linear shear and give rise to mean-field amplification. 

Another context in which the shear dynamo may be 
important is accretion disks, where turbulence is driven 



by the magnetorotational instability (MRI) Q . The MRI 
requires a (weak) large-scale field and gives rise to ve- 
locity and magnetic fluctuations at small scales. This 
turbulence could then conceivably couple to the large- 
scale Keplerian shear and amplify the large-scale field, 
thus closing the loop. While the feasibility of such an 
MPJ-dynamo mechanism will be the object of a separate 
study, we have verified already that the shear dynamo 
continues to work in the presence of rotation. 

We thank A. Brandenburg, L. Kitchatinov, G. Lesur, 
G. Ogilvie, J. Papaloizou, D. Shapovalov, D. Sokoloff and 
especially F. Rincon for valuable discussions. Simulations 
were done on UKAFF, NCSA (Illinois) and Cambridge 
HPCF. This work was supported by STFC (T.H., A.A.S), 
Newton Trust (T.A.Y, T.H.), US DOE CMPD (A.B.I.), 
the Royal Society (I.R., A.B.I.) and the Leverhulme Trust 
Network for Magnetized Plasma Turbulence. 



[2] 
[3] 
[4] 

[5] 
[6] 

[7] 

[8] 
[9] 

[10] 



[11] 
[12] 

[13] 



[14] 
[15] 
[16] 

[17] 

[18] 
[19] 
[20] 

[21] 



G. K. Batchelor, Proc. R. Soc. London A 201, 405 (1950) 
A. P. Kazantsev, Sov. Phys.-JETP 26, 1031 (1968) 
Ya. B. Zeldovich et al, J. Fluid Mech. 144, 1 (1984) 
M. Meneguzzi, U. Frisch, and A. Pouquet, Phys. Rev. 
Lett. 47, 1060 (1981). 

A. A. Schekochihin et al., Astrophys. J. 612, 276 (2004). 

A. A. Schekochihin et al., New J. Phys. 9, 300 (2007). 
K.-H. Radler and M. Rheinhardt, Geophys. Astrophys. 
Fluid Dyn. 101, 117 (2007). 

M. Ossendrijver, Astron. Astrophys. Rev. 11, 287 (2003). 
S. A. Balbus and J. F. Hawley, Rev. Mod. Phys. 70, 1 
(1998); J.-F. Donati et al., Nature 438, 466 (2005). 
K. T. Chyzy et al., Astron. Astrophys. 355, 128 (2000); 

B. M. Gaensler et al., Science 307, 1610 (2005). 

R. Monchaux et al., Phys. Rev. Lett. 98, 044502 (2007). 
I. Rogachevskii and N. Kleeorin, Phys. Rev. E 68, 036301 
(2003); Phys. Rev. E 70, 046310 (2004). 
E. T. Vishniac and A. Brandenburg, Astrophys. J. 475, 
263 (1997); N. A. Silant'ev, Astron. Astrophys. 364, 339 
(2000); M. R. E. Proctor, Mon. Not. R. Astron. Soc. 
382, L39 (2007); N. Kleeorin and I. Rogachevskii, Phys. 
Rev. E 77, 036307 (2008). 

E. G. Blackman, Astrophys. J. 496, L17 (1998). 

V. Urpin, Mon. Not. R. Astron. Soc. 308, 741 (1999); 

Phys. Rev. E 65, 026301 (2002). 

K.-H. Radler and R. Stepanov, Phys. Rev. E 73, 056311 
(2006); G. Riidiger and L. L. Kitchatiov, Astron. Nachr. 
327, 298 (2006). 

A. Brandenburg, Astrophys. J. 550, 824 (2001). 
E. N. Parker, Astrophys. J. 122, 293 (1955). 
G. I. Ogilvie, unpublished (1998); Y. Lithwick, Astro- 
phys. J. 670, 789 (2007). 

T. Elperin, N. Kleeorin, and I. Rogachevskii, Phys. Rev. 
E 68, 016311 (2003). 

A. Brandenburg, Astrophys. J. 625, 539 (2005). 

D. Shapovalov, private communication (2006). 

Y. Ponty et al., Phys. Rev. Lett. 94, 164502 (2005); P. D. 

Mininni et al., Astrophys. J. 626, 853 (2005). 

We found that the shear dynamo effect disappears when 

Re and Rm are too small, as predicted by SOCA [13J. 



